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III. 

On the Use of Equivalent Factors in the Method of Least Squares. 



By GEORGE P. BOND, A.M., 

ASSISTANT AT THE OBSERVATORY OF HARVARD COLLEGE. 



(Communicated, April 15, 1856.) 



One of the most important applications which has been made of mathematics to 
investigations in physical science has for its object to ascertain the best manner of 
combining data affected by unknown errors of observation, so that the probable effect 
of these errors shall be the least possible. The method of least squares proposes 
to accomplish this, by reducing to a minimum value the sum of the squares of the out- 
standing errors, and, by conforming to this single criterion, to fulfil the condition, so 
desirable in the prosecution of thorough and exact research, of reducing to its least 
possible amount the influence of errors in the data employed. 

The investigations here presented have been entered upon with the design of deter- 
mining the degree of numerical exactness proper to be observed in making use of the 
method of least squares, in order to secure its peculiar advantages with the least out- 
lay of labor. 

Some detail in the discussion seems to be called for from the prevalence of a practice, 
almost universal among computers, of adhering to the letter of the method of least 
squares with a strictness which implies a misapprehension of its true spirit. It is 
impossible to adduce any valid reasons to justify such a course when it must be fol- 
lowed at a serious expense of time and labor in the computations. 

It has not escaped the observation of Gauss, in his original exposition of the method, 
that some freedom of interpretation may be allowed when its theoretical results are 
applied in practice, as the following passage, referring to the solution of equations by 
least squares, will show : — 
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" When the number of functions or equations proposed for solution is considerable, 
the computations become laborious, the more so from the circumstance that the co- 
efficients by which the primitive equations are to be multiplied are almost always com- 
plicated decimal fractions. If it is not thought worth the trouble in such a case to 
calculate the products with exactness by means of logarithms, it will generally be suf- 
ficient to substitute for them (i. e. for the multiplying factors) more simple numbers 
differing but slightly from them." * 

In his subsequent researches, it does not appear that Gauss has given any further 
development to the suggestion here put forth. Indeed, the introduction of modifications 
of a like nature, however desirable in a practical point of view, would have deprived 
a purely theoretical discussion of much of its elegance and symmetry. Yet the passage 
above quoted lends the support of the highest authority to the leading proposition 
which we shall have occasion most frequently to insist upon ; namely, the propriety of 
allowing some relaxation of theory in applying the calculus of probabilities to the dis- 
cussion of data affected by ordinary errors of observation, whenever the modification 
conduces to convenience and the saving of labor at the sacrifice of no appreciable ad- 
vantages. 

Even an unqualified admission of the superior probability of results which exactly 
fulfil the criterion proposed in the method of least squares, does not relieve us from the 
necessity of restricting it to examples which never actually occur, that is, if the ques- 
tion be made a rigorous one; f — to such, for instance, as involve the discussion of obser- 
vations which are entirely free from unknown constant errors, or errors following any 
law of facility which does not imply the assumption that the mean error is proportional 
to the square root of the mean of the squares of the individual errors. But we know 
that this proposition, which lies at the foundation of the whole subject, is not suscepti- 
ble of absolute demonstration by any process of mathematical reasoning. Further than 
this, we know from constant experience that the law of distribution of errors recognized 
in the method of least squares practically fails, in extreme cases, both for very large and 
for very small errors. If any illustration of the failure of the assumed law be needed, 
it will be found in the familiar instance of computing by it the probable error of the 
arithmetical mean of a very large number of observations, where common sense assures 
us that the theoretical probable errors of the result are invariably smaller than they 
should be. 

Why, then, should an implicit adherence to its minutest details be required as essen- 

* Theoria Motus, § 185. t Theor. Comb., § 17. 
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tial to its successful application, or to the attainment of all the advantages which its 
employment may confer upon the discussion of any practical problem \ 

It is true that no other system can be proposed which is free from similar objections, 
or which can be mathematically demonstrated to be exclusively the best, without quali- 
fication, and therefore the arguments above stated are of no force whatever, if employed 
as reasons for the rejection of the method of least squares. They nevertheless greatly 
weaken the position of those who would insist upon a strict compliance with its pre- 
cepts, and effectually preclude all arguments of a purely theoretical character in support 
of such a course. Still it is desirable that the force of any objections which may be 
made to an attempt to modify the theoretical conditions for effecting the most favorable 
combination of equations should be appreciated at their true value. We therefore pro- 
pose to show that the spirit of the method of least squares, rightly apprehended, in 
reality rather invites than discountenances a liberal construction of its rules. 

Admitting that the best possible solution is attained when the sum of the squares 
of the outstanding errors, represented by il, is a minimum, it is evident that il is a 
minimum relatively to the manner in which the original equations have been treated. 
And since the peculiarity of the solution consists in the employment of a system of 
factors, a, a, &c, by which the original equations are multiplied before combination, 
the first differential coefficient of il relatively to either of these factors, in the case of 
the least-square solution, must have the value for each factor, 

dsi 

a a 

"When, therefore, the factors are varied by small amounts, 8 a, & a, &c., the conse- 
quent variations of il developed in a series, will contain only terms multiplied by the 
second and higher powers of 8 a ; or, in general terms, if we deviate from the exact 
(1.) precepts of the method of least squares by small variations of the first order, we shall fail 
to satisfy its fundamental criterion by small terms of the second order only. 

Looking thus at the most elementary principle of the method, we find a warrant 
for some degree of liberty in applying it, — a liberty which Ave can scarcely hesitate to 
avail ourselves of, if we further consider the peculiar circumstances attending its actual 
employment in the discussion of data furnished directly by observation. 

Among its first requirements is the assignment of weights to the original observations ; 
but it is one which it is not possible to fulfil correctly, for we are provided neither with 
a theory nor with data for the purpose. All that can be done is to accept, as indices of 
the relative value of the different observations, certain numbers depending either proxi- 
mately or remotely upon no other authority than the mere exercise of the judgment 
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alone. No one can pretend that this is a process susceptible of strict accuracy ; yet an 
error here is as fatal as if we had disregarded any other of the precepts of the method. 

This step being an arbitrary one, although one of fundamental importance, we may 
properly appeal to it as a precedent for the modification of others suggested by con- 
siderations of convenience, though they may not, like this, be justified on the plea of 
actual necessity. In this view of the subject, we find support for the modification 
suggested by Gauss, in the passage we have quoted above. Each of the complicated 
factors which it is there proposed to simplify is itself a product of two other factors, 
one of which is the weight of the equation under treatment; if one of these, that 
is, the number representing the weight, is erroneous, the product is of course errone- 
ous, with whatever accuracy the other is expressed. 

Again, as a matter of convenience, it is usual to express the conditional equations 
proposed for solution in a linear form, by reducing the indeterminates to small quanti- 
ties and neglecting the terms multiplied by their second and higher powers, and to con- 
struct from them normal equations, as they are called, previously to applying the 
method of least squares. Both of these may be practices perfectly allowable under the 
circumstances, but since they are almost always theoretically incorrect, their admis- 
sion is a virtual relinquishment of all pretensions to a rigorous course of computation, 
and cannot be compensated for by any subsequent refinements. 

We will now proceed to examine the limits of accuracy appropriate to the arithmeti- 
cal operations required in the combination of conditional equations by the method of 
least squares, and afterwards to develop in detail some proposed modifications of that 
method, haying for their object the reduction to its minimum value of the amount of 
labor requisite for its successful application. 

It is scarcely necessary to remark, that the subject is plainly one which is in its 
nature somewhat vague and insusceptible of rigorous treatment, though it is at the 
same time interesting from its practical bearings. If no very precise or definite rules 
for regulating the degree of numerical exactness suited to the discussion of any given 
problem can be arrived at, it may still be of service to point out the principles which 
ought to guide the computer in the choice of such limits as shall perfectly meet all 
reasonable requirements of accuracy, without imposing upon him the unprofitable labor 
of multiplying the extent and difficulties of calculation,. to no useful purpose, and 
without the remotest prospect of sensibly improving the real value of the results. 

Let us suppose a series of equations, 

a x -f- b y-j- -\- m = e, 

a x -f- b' y -f- -}- m = e ', 
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in which m is the element derived from observation, and e the unknown error of the 
equation, to be solved by the method of least squares, giving for x the value x, with its 
probable error, £, obtained from a comparison between the observed and the computed 
values of m, after substituting x, y, &c. in the primitive equations. 

If x be the true value of x, we may represent by x a quantity such that it is 
an even chance whether x Q — x is comprised between the limits e -\- x and e — x. 
The magnitude of the limit defined by x has an evident relation to the question how 
far the simplification of the arithmetical processes may be carried without detriment 
to the results. 

For instance, the solution of the above equations may be repeated with small varia- 
tions from the process at first applied, giving for x a new value x { with a probable 
error £, , differing but little from s. If we were in entire ignorance of the relative 
amount of the probable errors s l and s, there would be no reason at all for giving 
the preference to x rather than to x x . If only the single circumstance were known 
that Si exceeded s by a given small amount, we should be equally at a loss, while 
the value of s remained unknown, to state the relative weight of x x compared with 
x, and should, in fact, be again obliged to resort to the hypothesis that $ and « : were 
sensibly equal. And in general, the greater the uncertainty of s, or, in other words, 
the larger the value of x, the less reason would there be for excluding from competition 
with x any other determination of x, such as x u of which the probable error £j dif- 
fered but little from £. 

In order to employ the limit x as here proposed, its value must be known before the 
computations have reached an advanced stage. That this is not ordinarily practicable 
will readily appear. On the other hand, it must be left entirely to the judgment of the 
computer to decide as to the precise manner in which x is to be applied in limiting the 
allowable amount of difference ^ — «. 

Objections of a similar character apply equally to other standards which might be 
proposed for the same object. As has before been remarked, the question must be 
treated, if at all, upon a somewhat arbitrary basis, and we must be content with sug- 
gestions addressed to the judgment or common sense of the computer, in cases where 
no fixed rule is admissible. 

Viewed in this light, there will ordinarily be no difficulty in recognizing the point 
at which there will be danger of compromising accuracy in the attempt to simplify the 
computations, nearly enough at least for practical purposes, if we are prepared to 
admit, at least in its general spirit and tendency, the truth of the following propo- 
sition : — 
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The application of the method of least squares to the discussion of observations of physi- 
cal phenomena, with the exception of a few special cases of rare occurrence, requires the 
(2.) use of such numbers only, in the arithmetical processes peculiar to it and characteristic of 

the method, as may be designated by one of the numerals 0, 1, 2 9, or of the fractions 

— , — — , or by a product of one of these numbers by an integral power of 10. 

An idea may be formed of the amount of the intentional errors occasioned by these 
substitutions, by noticing that if by N is represented any number whatsoever, and by 
N' a number chosen from the proposed series which most nearly coincides with N, we 
shall have 

N— N' i 
/„ » The maximum value of — — = - nearly. 

, ,, , , N—N' i 
The probable value of — — <^ — . 

Before proceeding to a detailed investigation of the consequences of the changes pro- 
posed, it will be useful to point out the degree of insecurity attaching to the values 
which must ordinarily be adopted to represent the probable error of x ; the different 
sources which may be supposed to contribute to the increase of £ ; and their relative 
importance in connection with the question of the comparative accuracy of the two 
results x and x l . 

a may be referred to the combined influence of two mutually independent errors 
7) and t) , t) being the probable value of x — x which would result from the errors of 
observation alone, supposing the theory of the method of least squares and its applica- 
tion to the data to be rigorously exact, and rf the probable amount of error in x having 
its origin in errors necessarily committed in the discussion of the observed data, sup- 
posing the mode of discussion, although the best practicable, to fall short of strict 
conformity with the theory. r)[ represents the value of r/ when the same data have 
been reduced, by a process made intentionally still less exact, to a small extent, both 
in its theory and in its arithmetic, than that which gives the error r/. rj[ will bear 
to x L a relation similar to that which rj bears to x. y cannot be completely eliminated, 
so long as the errors of observation remain unknown, by any treatment, and the same 
may be said of if ; but t)[ can always be reduced to its least limit, if, by suitable refine- 
ments of theory and of computation. In view of the fact that i) and rf must have 
always sensible, but very uncertain values, it will be of but little consequence that 
r][ should be reduced to its utmost limit without regard to the labor and inconvenience 
which it may cost. At all events, the attempt will be ineffectual as a means of im- 
proving the substantial accuracy of the results, as we shall presently see. 
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Since i\ is independent of if and */, > we have, assuming if to be the least attainable 
value of %, 

If if' be used to designate the probable value of 00 — x l which would result from 
small intentional deviations from that treatment of the data which is recognized to be 
the best, we have 

(4.) ,',■ = ,' 9 + V*. ? = r? + r!\ „« = ,2 + ,'2 + ,»». 

As regards the uncertainty of £, some estimate of its extent may be obtained in the 
following manner. 

If it is an even chance that the error of which the probable value is i) is comprised 
somewhere between the limits »? -j- X and i) — X, 17 having been derived from compari- 
sons of a given system of equations with observation, the number of individual equa- 
tions thus compared being represented by n, and the number of unknown quantities 
entering into them by ri, X may be found from the expression* 

(5.) x = 0.477 



\/n- 



Any value of n — ri less than 100 gives 



-^ 21 ' 

The scale of substituted numbers (2) admits, as we have before stated, of represent- 
ing 17 within the probable amount of — y; hence, for any value of n — ri less than 100, 
the series will afford numbers representing i\ with a probable error less than X. A 
slight examination will show that a similar remark applies still more decisively to s. 

The considerations which oblige us to attribute a sensible value to if are too many 
and too obvious to require to be specified in detail. It will be sufficient to cite one or 
two which have already been alluded to. The existence of unknown constant errors in 
the data will render the application of the method of least squares, strictly speaking, 
inexact. From this source if will inevitably acquire some influence. Again, the un- 
certainty incident to any attempt to assign to the original data their proper relative 
weights, will have a similar effect. No process more loose and arbitrary can well be 
conceived, than that by which the relative precision of the elements afforded directly by 
observation is graduated. Yet, imperfect as it is, improvement in this particular is 
scarcely to be hoped for. Exact conformity with a theory which requires a previous 
knowledge of the relative weight of observations is quite impossible. 

• Gauss, Zeitschrift fur Astr., B. I. Theor. Comb.,§ 40. 
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At the same time, then, that the existence and influence of if are admitted, its amount 
is altogether uncertain, to an extent sufficient at least to make the uncertainty of s 
which is dependent on that of ij and if not less in proportional amount than that of i\ ; 
consequently we shall obtain from (5) the expression 

(5 a.) x > 0.477 



by which to measure the uncertainty of e. If n — n' < 100, 
(5 b.) * > - s. 

When t is represented by a number chosen from the series (2), the probable error of 
the representation is, by (3), less than — s; in other words, it is more than an even 
chance that this number will fall within the limits s -f- — s and s — — s ; and since 
the inherent uncertainty of e makes it more than an even chance that its actual value is 
outside of the limits a -j- — e and s — —s, in accordance with the above determination 
of x, we conclude that s can be represented by one of the series of numbers 0, 1, 2 

..... 9, or of the fractions -, - -, or by a product of one of these numbers by 

an integral power of 10, with more accuracy than we can determine its amount by 
one hundred comparisons between the observed and the computed values of m. It 
would be easy to show, from the probable existence of constant errors alone, that an 
indefinite increase of the number of comparisons with observation would not sensibly 
diminish the uncertainty of s below the amount stated. The proposition (2) would thus 
be sustained, as far as relates to all expressions for probable errors and weights, since 
they must depend upon conditions similar to those limiting the accuracy of s. 

An immediate consequence of this admission will be the extension of the proposition 
in question, in the qualified sense, at least, in which alone it is to be understood, to all 
other arithmetical expressions required in the application of the method of least squares, 
since the peculiar province of the latter is restricted entirely to the solution of equations 
of the form 
(5 c.) a (x — x t ) -f- b (y — y,) -f- + (m — %) = e, 

in which each separate term and factor may be defined as proposed in (2). 

To illustrate this, let us suppose for the moment that x l has been derived from the 
same primitive equations, but by an essentially different process from that by which 
x has been obtained ; x l would still be precisely equal to x, if it were not for the errors 
e, e', &c. Any such process, not intentionally bad, must evidently lead to a determina- 
tion of x x differing from x by an amount of an order not higher than that of s, while 
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one adopted expressly on account of its good qualities, though without bearing any in- 
tended resemblance to the method of least squares in its characteristic features, will 
diminish the difference x — x x to a value either less in absolute amount than s , or very 
nearly equal to it; so that (x — x x ), (y — y x ), Sec. may be sufficiently well expressed 
by the series (2) in the same sense that « may be. The equations 

ax -\-iy -\- -f- m = e, 

a x x + h y t -f- + mi = 0, 

in which e is the value of e obtained by substituting x, y, Sec. in the primitive equa- 
tions, give at once 
(5 c.) a (x — »,) + b (y — y t ) -f + (m — m,) = e, 

The solution of which by the method of least squares gives (x — x t ), (y — y-^) , 

and thence 

X = #1 -[- (x «i), 

y = 2/1 + (y — sa), 

Therefore the application of that method may be confined to the equations (5 c) alone. 

But by what has already been said, (x — #,), (y — y x ) , and e, are quantities 

sufficiently well expressed by the series (2), from which we readily infer that there can 
be no appreciable advantage in giving to (m — m { ) or to the products a(x — x x ), 
b(y — yj, Sec, or to the factors a, b, Sec, any higher exactness of expression. Hence 
the proposition (2) may be extended to every term and factor of the equations (5 c), 
and therefore to all the numerical processes peculiar to, and characteristic of, the meth- 
od of least squares. 

This extreme application cannot, however, be recommended even on the ground of 
convenience or simplicity ; on the contrary, the indiscriminate use of the fractional 
terms of the series would often be highly inconvenient ; and a form of solution like 
that just indicated would not always be desirable. 

In considering the different sources of error from which e and s x acquire their value, 
it will be convenient to compare the increase given to s by the introduction of small 
intentional inaccuracies, in consequence of which s becomes s lt either with s itself or 
with s — rj, by means of the following relations derived from (4) : — 
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Used in connection with (3) and the limit defining the uncertainty of s, 

E 1 

x > 0.477 . , or usually x _> — t. 

V» — n' 21 



Si will thus be comprised between the limits s -f- k x and s — k x, when if' has such a 
value that 

\e J 10 n — ri ' 

The relative accuracy of x and x l will now be investigated for some special exam- 
ples of deviation from a strict compliance with the method of least squares. 
Let the equations proposed for solution be the following : — 

ax-j-Jy-j- -j- m = e, weight = w, 

(7.) a 1 x' -f- b' y -f + ml = e', " = «/, 

where e is the difference between the observed and computed value ofm; m being the 
element derived from observation. 

In solving these by least squares, the final equation for x is formed by taking the 
sum of all the equations after multiplying the first by a iv, the second by a! w', and 
so on, and then making 

a w e -J- a 1 w' e' -j- =0, 

and fory 

b w e -J- V w 1 e' -)- =0, 

continuing in succession to form new equations until a final equation is obtained for 
each unknown quantity. 

We shall compare the results of two solutions of the above equations (7), in one of 

which (I.) the factors aw, a' w' conform strictly to the method of least squares. In 

the other (II.), these factors are replaced respectively by a, a! ; a being that one of 

the numbers 0, 1, 2 9, or of the fractions —, — — , or of their products by an in- 
tegral power of 10, which approaches most nearly to a given ratio with a w, and a that 
which approaches most nearly to the same ratio with a' w', fyc. In a similar manner, 
/3, 7 are used in the place ofbw, civ 

The true values of x and y we will indicate by x , y Those deduced by 

(I.) will be denoted by oc, y , and those deduced by (II.) will be denoted by #„ 

yi For the final equation for x, we make 

aw e -\- a' to' e' -\- =0. 
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For the final equation for oc x 

«e 1 + n'« 1 ' + = 0. 

For the corresponding final equations for x , which must be rigorous, we make either 

a w e -\- a' w' e' -f- = «tnc l -|-fl'i(i'e' 1 -f 

or 

« e + »' e 'o + = ae -\- a' e' -{- 



according as the first (I.) or the second (II.) form of combination is adopted. 

e , e, and e t are the values of e when the indeterminates oc , y Q , oc, y 

x x ,y x , &c, replace x, y in (7). 

The final equations for the combination (I.) are : — 

P x + P' y + P" 2 + + L = 0, 

(8) P l x+Qy-\-Q>z + + .M=0, 

{ '' P"x + Q'y + Q"z + + N = 0, 



P = w a a -j- w' a' a' -\- , Q = wbb -\- w' b' b' -f- , R = w c c -{- w' c' c' -\- 

(9.) P = wab + w'a'b' + , Q' = w b c + w< V c> + , R 1 = w c d + w'c'd' -f 



From the conditions I. and II. applied to the original equations (7), if we make 



Aa=aw-\-da, B$=bw-\-h$, Cy — c w -\- by , . . , 

A a' = a'w' + ba', P/3' = S' K>' + 8 /3' , C y 1 = c'w' + by 1 , . . . 



may be obtained 

P(tf- a;i ) + -P(3/-2/0 + + L + Px l + Piy l + = 

£-t-Px 1 + P'y' 1 + 

Hence, 

P (* — *,) + P (y — y.) + = 8« e, + ba' e\ + . 

(10.) P( x - Xt ) + Q (y_ yi ) + = S/3 ei -j-8/3'e' 1 + 



And in a similar manner, 

P (*„ — x) + P (y — y) -f = a Mi e + a' w> e' + 

(11.) p>( Xo — x ) + Q ( yo _ y ) + = Jtoe +J'w'e'o+. 



Since s is the probable value of oc — #, and ??" the probable value of {oc — oc,), to ob- 
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n 

tain the ratio — , we will compare the probable values of x — x and x — x u having, as 
above, 

/ j„ \ v" Probable value of (a; — »i) 

* Probable value of (x — x ) 

x — x and x — x i must be derived from a solution of equations (10) and (11), but 
since (II.) differs from (I.) by small variations only, we have, very nearly, 



(13.) w e, 2 + w 1 e', 2 + = w e s + w' e' 2 + . 



For the second member of (13) is a minimum relatively to the mode of solution, and, 
as has already been shown, (1), it differs from the first member by small terms of the 
second order only, those of the first order vanishing with the first differential coefficient 

of n = we 2 + w'e' 2 + 

If, then, fi , ix, and /^ represent the probable values of e , e, and e x corresponding 
to the unit of weight of the equations (7), we may assume, for the purpose of deter- 
mining x — Xi, that /J. — /itj is a small quantity compared with fi, since we have 

P w e 2 4- w 1 e' 2 4- 

~ = — an — TTa^i = *> verv near 'y- 

Pi vi e{ -\- w 1 e'f -j- 

Moreover, in the absence of exact knowledge of the magnitude of the errors of e , 

e' , , it is necessary to admit that they are best represented by the errors e, e', ; 

hence we have — = 1, and consequently - = 1, very nearly. 

The conditions of the solution (II.) give for the probable value of either of the ratios 

8_« 8_/3 8 a' 8 /3' 
aw' bw'' ' a'w' ' b'm 1 



(14.) >JL = 8 , {£ = ,, *lL- 8 , lE- g , 

n. VI h m n' »#»' hi tnl 



Because a being by (II.) nearly proportional to a w, the probable value of B a will also 
be proportional to aw ; and a similar remark applies equally to B 0, B y, &c. 

The probable values of the second numbers of (10) and (11) are then, respectively, 

8 a e x -\- 8 a' e'l + = g Pi V? » a w e ~\- a' w' e' -\- = ^ ^/p ? 

80 e, 4- 8/3' e', 4- = g ^ </q, b w e + b> w> e'„ + = Mo y-Q, 



Hence, in consequence of the identity of the coefficients P, J" , Q, Q' , &c. in 

the two systems (10) and (11), we obtain 
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Probable value of (a; — x { ) ix. Probable value of (y — yA u^ 

(15.) - = _ g — g; — — = — g = g 

Probable value of (a? — x) fi ' Probable value of (y — y) /*o 

And from (12) the general expression 



(16.) 


1 = * 

s 


Giving to g the value (3) 








we shall have in the present case 






V' i 




7 < ii' 


and by (6), 






^ 1250 



In other words, by using the form of solution (II.) in the place of a rigorous appli- 
cation of the method of least squares, the probable errors of the concluded results will 
not be increased by one one-thousandth part, — a difference entirely too small to be 
sensible. The two processes, as far as regards accuracy, therefore, may be consid- 
ered as perfectly identical. On the other hand, the advantages of simplicity and 
convenience are altogether in favor of the second, in which all the operations of mul- 
tiplication and division required in the construction of the final equations are reduced 
to their simplest arithmetical forms. 

The necessity of distinguishing between the probable error of <a? IS that is, the proba- 
ble value of (x — xj, and the difference between x x and x, or (ac l — x), deserves par- 
ticular attention here. "While x x will often differ very much from w, this fact taken by 
itself by no means indicates that the chances that x { is the true value are not sensibly 
as good as that x is. The discrepancy really proves that the original observations upon 
which the discussions have been based are so faulty, that very little confidence can be 
placed in either result, or in any other that can be deduced from the same data. 

To give completeness to the investigation, we will compare the processes by which 
the probable errors of the values of the indeterminates in the two solutions (I.) and (II.) 
are obtained. 

For the system (I.), let x, y, s, Sec. be eliminated in succession from the equations 

(8), in the following manner. 

P' 

Multiply the final equation for x by — and subtract it from the final equation for y, 

P" 
and again by — and subtract it from the final equation for s, forming the new equa- 
tions, in which x does not enter : — 
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(Q-^P')y+(Q'-|P")*+ + M-|'L=0, 

f 17 1 P" P" P" 

1 ; (Q!-j-P')y + (R—pP")* + + n-^l = o, 



to which we shall give the following notation : 

Q,y+Q',% + + JMT. = 0, 

(18.) Q' x y + R x z + + JV I= 0, 



Let# be eliminated from these new equations (18) in the same manner that x was 
from the equations (8), multiplying the first by -^ and subtracting it from the second, 
giving 

(B. - ?? Q>.) z + + N,-2jM, = 0, 

with the corresponding notation, 

(19.) *.,* + + N t> = 0. 

For a more complete illustration of the notation, we will collect in one view the 
equations obtained by the successive eliminations : 



(a) P x + P'y + P" z + + L = 

(b) P x+Qy+Q' z + + M = 

(c) P" x -\- Q' y -\- R ?-f + 2V=0 



Final 
Equations. 



(20.) (jb.) Q x y+Q', Z + + iM;=0] . 

( C .) Q^ + BxS + + tf. = lEquahonsformedby 

eliminating x. 

(c xy ) Rxyt -\- + N zll = i Equations formed by 

. J eliminating x and y. 

Let ft, be the probable error of one of the original equations of the unit of weight. 
The probable errors of the second members of the equations (20) will be, 

Probable error of the equation (a) = /x ^/p 

(b) =mVQ 

(c) =^vs 



(21.) " " " (*,) =fVQ. 

(O = /*y% 
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When the final equations obtained by (II.) are solved in an analogous manner, and 

the notation is changed so as to indicate the coefficients of x x , y l5 &c, P n P' l , 

which replace P, P\ &c, we have 

(a) P, x, + P, y, + P", *,+ + i, =0-J 

(i) 'Q.*.+ Qi yi+Q', *,+ + Jf, =o[ Final 

(c) "R l ^ + 'E, y, + B, «, + + JVi = j Equations. 



(22.) 



(23.) 



(».) 
(«.) 

(W 



Qi.yi+Q'i. «i + + J«r„ =o] 



'Bi.yi + ili. *i + + ^i. =0 



Equations formed by 
eliminating «,. 



Bi ty Z\ -|- — j— -ZV ljE! , — •» Equations formed by 

) eliminating x L and y,. 



Probable error of the equation (22) (a) = (1 ± - g 



(») =(!±i« 



(0 



(i ± s * 



(».) =0±i<r 



(c.) = (l ± J *; 



\r, 



for. 

"sir* 

IB,,; 



<o =( 1 ±i^)/'JI^' 



To demonstrate these results, (23), it is to be observed that the probable errors of the 

second members of the equations (22), (a), (b x ), are the probable sums of the 

second members of the equations 

P. fa — *i) + -P'i (y« — yi) + = « e„ + «' e'„ + 

Pi 



Qu (y. — yi) + = — ^ «) e + (/* - '■%■ a') e'„ + 



Pi 



The probable value of a 2 eo is, by (13) and (14), 



.4 a 



2 2 f i i & a \ tt a a /i . x a a 
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and the probable value of the sum of the terms a e -\- a e' -\- is 

« e„ + «' e'„ + = -J| V(TT^) V«« + «'«' + = m (1 ± I g) J J • 

which is the probable error of the equation (22), (a). 
Again, since 

Aa = au,(l + —\, B/S = i»fl + ^.V 
\ awj \ bwj 

and 

P t = « a 4- «' a' -f , Q t = p b + 0' 5' + 

P'! = a b + a' J' + , 'Q t == /3 a + /3' a' + 

or, substituting the probable values 
we have 

(,_a.y < ., (I±I)[ i(„_s.,) + .(3.._^. )] . 

Moreover, 

The sum of all the terms ('S& a a — ^ a\ = f^? P. — ^ '©A = 0, 

\pi Pi J \ii Pi J ' 

Therefore the probable sum of the terms 

(0- : |«)eo + (/3'-^«')e'o + 

X 1 Xi 



will be 



vnnJf=-..<i±J,>Jf--. 



which is the probable error of the equation (22), (h x ). 

The other probable errors in (23) are readily supplied by analogy. 

If we neglect - g, of which the probable value is less than — , the probable errors of 

2 50 

(22), (a), (b x ) become 

Probable error of equation (22) (a) = ^ I— , Probable error of equation (22) (5 r ) = p I— . 

We shall now proceed to explain a third form of solution, (III.). 
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Returning to the equations (10), 

P (* — *,) + P' (y — y,) + = «i»« +e' 1 8a' + , 

(24.) p/ (a ._ a . i) + Q { y — yi ) + = ei8 ^ +e / i8/3 , + i 



we find, for the probable values of their second numbers, 



g /x Vw a a -\- w' a! a! -\- = g ^ \/P , 

(25.) g /u Vwi b + w' b' V + = g p, VQ, 



It is evident that the probable sum e x 8 a -f e\ 8 a -f- , being proportional to the 

square root of the sum of the squares of the individual terms, depends mainly upon the 
large terms ; or, since e L 8 a = e x Vw— = and e { vV= p, this sum will be 



Vio 



la« 8 »«■ 

N io w' 



r + 



If any two or more of the coefficients -j=, as, for instance, -j=- and 77%, were equal, 
any small change increasing the former and diminishing the latter by equal amounts 
would not alter the coefficient of /x; but if -7=- were much smaller than -7^-, we 
should have, very nearly, 

So 8 8«^ _ So 2 

V) ' ID W ' 



and a small change in r^=~ would affect the coefficient of p by an amount insensible 
compared with the effect of an equal change in —j= . 
Let (P) represent the sum of a certain number of the largest of the terms composing 

the series 

P = w a a -\- w' a 1 a 1 -f- 



and (p) the sum of a number of the smallest of the terms of the same series. Let also 
(8 P) be the sum of the terms — corresponding to the series (P), and (8p) the sum of 

the terms — , corresponding to the series (p). 

Then we have the probable values 

For the large terms, (8 P) = g* (P), 
* "' For the small terms, (dp) = g* (p), 

g representing the general probable value of — for all the terms, whether of large, 
small, or medium value. 
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Let us suppose the mode of solution (II.) to be itself varied by changing the factors 
a, a', &c., corresponding to the large and small terms, so that for the large terms, g, or 

the probable value of — , -?-,, , for these particular terms, becomes 

■"■ aw a w '■ 

and for the small terms, 

g = h. 

We shall then have the probable values, 

For the large terms (8 P) = IF (P), 
For the small terms (Spi) = W (p). 

(8"P) and (8p) becoming (8 P,) and (8pJ when g becomes H and h. 
In order that the probable sum of the second member of the equation, 



P (c — *,) + P {y — Vl ) + = e, 8 a + e\ 8 «' + 

should not be increased by the proposed changes of S a, we must have 

or, by (26) and (27), 

fl 2 (P) + A 2 (p)<g 2 (P)+^(p). 

We shall assume, for the terms corresponding to (p), that the probable value of h is 

h = — 1. 

This condition involves only small changes in the factors a, a , because, for the 

terms corresponding to (p), a w being small, 8 a = aw h = — aw, will also be small ; 

we then have 

H°(P)<g>(P) + (g>-l)( P ), 

or, since we can put g % — 1 = — 1 very nearly, g being small compared with unity, 
we obtain 

(28.) g *_iP>M, iP<g*-M, 

(P) (P) 

representing the condition to be observed in order that the second members of (24) 
should not be increased by the changes made in the large and small values of «. 

This, it will be remembered, can be applied only when the condition h = — 1 in- 
volves only small changes in the factors a, a' of the order of the mean value of 

8 a for all the factors. yQ being necessarily a positive quantity, H must always be 
less than g. 
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(28) may easily be extended to the analogous cases of the second members of the 
equations 

P' (x — *,) + Q (y — y.) + = e 1 8/3 + e' l 8/3' + 

P" (» — xi) + Q' (y — yi) + = My + e',8/ + 



so that 

(29.) g»_2F>ilL, ?-&>" 

give the limits within which the proposed changes of the factors /3, /3' 7, 7' 

will not increase the probable sums e x 8 /3 -|~ e\ 8 ft -J- and e x 8 7 -(- e\ 8 y -f- 

For the factors a, a' , @, $ corresponding to the equations most important 

in their influence upon the final determination of #, y respectively, if we use num- 

bers chosen from a series for which — — — is only — as large as it is for the series (2), 



N—N 
ova. a series tor wmcn 

we shall have 



10 8 



And if at the same time we omit altogether a certain number of the unfavorable 

equations by making in these instances a = 0, /8 = , that is, S a = — aw, 

S/9 = — b w, or h = — 1, we find 

5 631 

We shall therefore keep within the limits (28) and (29) as long as the coefiicients in 
the omitted equations satisfy the conditions 

M<-_L (g) < 1 

(P)*^63l' (Qj^esi' 

The probable values of x — x^y — y x , will not have been increased, and conse- 
quently the solution may be accepted as equivalent to II. 

A general method, III., of adjusting the degree of numerical accuracy which should 

be observed in the expression of the factors «, a' , 0, ft , may be derived 

from the following considerations. 

In II. the adjustment is evidently not so favorable as it might be, since the limit of 
the intentional inaccuracies 8 a, 8 a' , 80, 8ft has been fixed by the relations 



8 a = a w g, 8 a' = a' w 1 g', dp = bwg, 8 £' = V w' g 1 , 



g having the same average value whether aw, bw be large or small ; thus the 

VOL. VI. NEW SEEIES. 26 
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largest inaccuracies are committed when aw, bw are largest, that is, when the 

equation has most influence upon the final result for any particular indeterminate. 

In order to secure a more advantageous distribution, it will be necessary to recur 
again to the equations (10). It appears that, for a given limit of inaccuracy in the 

expression of the factors, the probable values of x — x x , y — y x will be least 

when the separate terms of the second members of these equations, irrespective of their 
signs, are equal to each other, or 

e t 8 a = e'l 8 a' — e x 8 /3 = e\ 8 /3' = 

8 a, 8 ft ought therefore to be inversely proportional to the probable errors of the 

primitive equations, or directly as the square roots of their weights. 
We shall, then, define III. by the relations 

8 a = A g Vw~, 8/3 = Bg Vw, 

(30.) 8a' = A.gVw', &p'=BgVw', 



A, B and g being constant quantities. 

To secure the equivalency of II. and III., the values of A, B must depend on 

the condition that the probable values of the second members of (10) should remain 
unchanged, or 

S,*VA 2 + A S + = g?\/P, g M VB 2 + B 2 + = g^VQ, 

AVn= Vp, bV»= Vq, 



Hence it is easy to conclude, that, if we make in (30) 

A = mean value of a \/w, B = mean value of b Vw, 



the means being in every instance taken without regard to signs, the probable values of 

x — <#i, y — y x will be smaller in III. than in II., while III. in point of facility 

has a decided advantage over II. ; since by making a = 0, ft = in all cases in 

which a \/w < A g, b i/w < B g, , a considerable amount of unnecessary computa- 
tion may often be avoided. 

The following will then be the limits of intentional numerical inaccuracies allowed 
in the expression of the factors a, a' , /S, /3' in the three methods. 

I. II. III. 

8a =0, 8/3=0 da=awg,8p=bwg 8«= kg Vw, 8/3 —Jig Vw 

8 u' = 0, 8 £' = 8 a' = a' w' g, 8/3' = l' w 1 g 8 «' = A g Vw', 8 /3' = B g \/w' 
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The following tests of the correctness of all the numbers entering into the final equa- 
tions will be found useful. They apply equally to the three methods. 
Let the sums S, S' , 2, 2' be formed as follows : — 

«+&+ + m = S , «+/ 3 +y+ = 2, 

( 81 -) a' + J'+ + m' = S / , «'+j3' + y+ = 2>, 



We shall then have the following test controlling at once all the computed quan- 
tities, P 15 P\ , 'Q 19 Q x , L lt My 

2 (P» Qi, Bi ii, JUi, M ) = J? (Si). 

For we have 

P, + P'i + + L l =« a +«5+ 

+ «' a' + a' J' + ^ :^ S a + S' a' -\- 

'Qi + Qi + + Mi = /3 a + 13 J + 

(32.) + ff a' + /3' 5' + }> =SHS'^' + < 

+ 

"Bi + 'Ri+ + JV, = y a+ y J+ 

+ ya' + y'5' + }> ^Sy + S'y' + 

The sum of all the equations (32) is 

(Py + P'^ + £,) + (<Q, + Q, + + Jf,) + = 8^+8'^+ 

or 

(33.) 2(Pi,Qi,Rt £» a*i, JVi ) = J£ (Si). 

When this is not fulfilled, the particular equation at fault may be detected by using 

the partial tests, 

Pi + P'. + In = S « + S' «' + 

(34,) 'Qi+Qi+ M 1= S0 + S'/3' + 



If all the numbers are not carried out to the same number of figures beyond the 
decimal point, it will be advisable, in order to apply the test to the best effect, to alter 
arbitrarily the decimal point, making, for the time being, the number of figures to the 
right hand of it the same in all the equations. 

The successive stages in the computations by (II.) and (III.) will be as follows : — 

a) To assign the values of the equivalent factors a, ft , which can readily be 

done by simple inspection. The factors once adopted must be used rigorously in mul- 
tiplying every term of the equations to which they belong. 
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h) After the multiplications have been performed, and the sums taken, the numbers 
adopted in the final equations are to be tested by (33). 

c) The solution of the final equations. 

d) The determination of weights. 

If changes have been made in the decimal pointing, or otherwise, by introducing 

the constants A, B , it must be remembered that, although the final equations 

thus formed will give the same values of x u y x , &c. that would have been obtained 
if no such alteration had been made, the determination of the weights and probable 
errors of x x , y x , &c. requires that the correct pointing be restored in the coefficients, 
or else that the probable errors be computed in conformity with the formulae (23). 

When the number of indeterminates is considerable, it will be advisable, in solving 
the final equations, to eliminate x, y, z, &c, in succession, and then to repeat the 
operation, commencing the elimination in the reverse order, z, y, x, &c. One of the 
advantages of so doing is a complete check upon the work by the comparison of the 
value of that indeterminate which is obtained last by both eliminations. It is, however, 
mostly recommended from its facilitating the computation of the weights. In this case, 
the following formulae may be used, if the number of indeterminates does not exceed 
six. Let these be x, y, z, !-, y, f, and their weights, W (x) , W (!/) , &c. The ordinary 
formulae for computing the weights give 

(35.) PF(0 = TJ X y x £ , , W ix> = P ? , { , y . 

W(q is the coefficient of f in the equation resulting from the elimination of x, y, s, f, 
and 17, by the process indicated in (18), (19), and (20), and W ix) the coefficient of x 
after £, rj, f, z, and y have been eliminated. We have, also, 



(36.) W w = %JLU W (0 , W w = fcU! JTc 

u x y z l ^; <t I z 



The factors and divisors required in (36) will have been already computed during 
the eliminations which have preceded. 

From the equations containing only f, 17, f, the latter is to be eliminated ; and from 
the equations containing z, y, and x, x is to be eliminated. We then have 

^* y z ? tit wr &C « l x 



Wrs — y " ? Wr ^ 

J-x y z C 



W^ = ^^W ( yy 



y z (, %«{* 

For the weight of a? l5 y { , &c, when (II.) or (III.) is used, we shall have, from (23), 

Weight of «, = AW (xy , 
(37.) « Vl = B W w , 
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The limits of effective accuracy appropriate to the numerical representation of the 

coefficients a, b may be investigated in the following manner : — 

If we determine x x , y i by the method of least squares from the equations, 

(38.) (a — d a) x t -|- ( b — d b) y x -)-. -\- m — d m = e x , weight = w, 



we shall have from (8) and (38), if d a, d b , which may be employed to represent 

the errors of the adopted coefficients, are small, 

(39.) P ( Xl ~ X )+ P '(^ — ?) + = aw {— e-{-xda-\-ydb + + dm J + 



If // is the probable value of dm Vw, the most suitable values of x d a, y d b . . 
evidently fulfil the conditions 

(40.) xda\/w=ydb*/w = dm\Sw = fi', xda' Vio' = ydb'Vw' = dm'\ , w' =/*',. 



Observing that we may substitute in the second members of (39) the probable values 

e < a x, e' < a' x , e < by, e' < V y , we may conclude, by comparing (39) 

with (11), (12), and (6), that, if/*' is less than the limit 

(41.) s '» 



/*' = 



yV + 2 

ri denoting the number of indeterminates in (38), the difference s t — s of the probable 
errors of x u y x derived from (38) and of x, y derived from (8) will be 

(42.) , 1 _,<ig/» i . 



2 



Consequently, if g < — , £, — s will be less than — s or less than in II. 

When p is known, the limits for admissible values of da, db in (38) will be 



I* U, IX! X 

(43.) da<C — , da' <^ — Kl — = db = -da, db 1 = -da' dm = xda, dm' = xda' 

xVw a;vv y y 



If the mean values of ax, by , irrespective of their signs, are all of the same 

order of magnitude, we may substitute in (43) the a priori probable values 

(44.) x<f rz» , ^=? 

A \n + n' — 1 y A 

where n is the number of equations (38), n' the number of indeterminates, and, taking 
the means without regard to signs, 

(45.) A = mean value of a \/w , B = mean value of b \/w M = mean value of m VV. 
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We will now proceed to compare the three methods of solution (I.), (II.), (HI.), by 
applying them to the following series of equations of six indeterminates, taken from a 
memoir, by Gauss, on the elliptic elements of the orbit of the planet Pallas.* 

Original Equations.^ 

1 = — 183J3 + 0.79363 d L + 143.66 d y + 0.39493 d w + 0.95920 d <p — 0. 18856 d Q + 0. 17387 d i 

2 = — 6.81 — 0.02658 di + 46.71 dy + 0.02658 djr — 0.20858 d v + 0.15946 da-\- 1.25782 di 

3 = — 0.06 + 0.58880 di + 358. 12 dy + 0.26208 d ■n— 0.85234 d 9 + 0.149 12 da + 0.17775 di 

4 = — 3.09 + 0.01318 di + 28.39 dy — 0.01318 dit — 0.07861 d <p + 0.91704 d a + 0.54365 di 

5 = — 0.02+ 1.73436 di+ 1846.17 dy — 0.54603 d*r — 2.05662 dg> — 0.18833 da — 0.17445 di 

6 = — 8.98 — 0.12606 di— 227.42 dy + 0. 12606 dn- — 0.38939 d 9. + 0.17176 do — 1.35441di 

7 = — 2.31 + 0.99584 di+ 1579.03 dy + 0.06456 d,r+ 1.99545 d 9 — 0.06040 dQ — 0.33750 di 

8 = + 2.47 — 0.08089 d L — 67.22 d y + 0.08089 d * — 0.09970 d 9, — 0.46359 d a + 1 .22803 d i 

9 = + 0.01 + 0.65311 dL+ 1329.09 d y + 0.38994d,r — 0.08439 d 9 — 0.04305 d a + 0.34268 d t 
9, = + 38.12 — 0.00218 di+ 38.47 dy + 0.00218 d„- — 0.18710 d q, + 0.47301 dQ — 1.14371 di 

10 = — 317.73 + 0.69957d L+ 1719.32 dy + 0.12913d7r — 1.38787 d9> + 0.17130dQ — 0.08360 di 

11 = + 117.97 — 0.01315 di— 43.84dy + 0.01315 d,r + 0.02929d 9,+ 1.02138 da — 0.27187di 

The probable error of one of these equations is /* = + 90", the weights being 
equal, excepting for 9 a , which has been excluded from each of the solutions. 

As an illustration of the mode of applying the limits (43), we will make in (41) and 
(42) g' = — , we shall then have 

M = ± 90", n = 11, A=±0.5, M=±60", (x = d i)< ± 100" , da<± 0.013, 
g'= ^, n'= 6, B=±700, /=±1".8, (y = d y )< ± 0".07, dJ<±18. 

0, D, E, and F, being of the same order of magnitude with A, we may conclude from 
(42), (43), and the above values of d a and d b, that, if we reject in all the equations the 
two right-hand figures from the values of m, and the three right-hand figures from all 
the other numbers, writing, for instance, the first equation 

= — 184" + 0.79 d i + 140 d y + 0.39 d n + 0.96 d <p — 0.19 d a + 0.17 d i, 
and the second 

= — 7" — 0.03 d i + 50 d y + 0.03 d n — 0.21 d <j> + 0.16 d a + 1.26 d i, 

the probable error of the value of either of the quantities of dL, dy , when the 

equations thus written are solved by the method of least squares, will exceed the prob- 

* Disquisitio de Elementis Ellipticis Palladis. Comment. Soc. Reg. Gottingensis Eecent. Vol. I. 
t The computations necessary for the solution of these equations have been executed by Messrs. J. F. 
Flagg and T. H. Safford, jr. 
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able error of the same quantity obtained with the employment of the exact coefficients 
by less than — of that error. 

J 1250 

We shall, however, for the present confine our attention to a direct comparison be- 
tween the results of the solutions I., II., and III., and retain in each the exact coeffi- 
cients of the original equations, adopting the constants A, C, D, E, F = 1.0, and 

B = 1000.0, and for the factors aft a somewhat ruder system of representation, 

especially in II., than that upon which the previous discussions have been based. 

If we employ the limits (28), (29), we find that the most favorable equations are (1), 
(5), (7), and (10), and the least favorable (2), (4), (8), (11). Applying (III.) in the 
particular form to which it is limited in (28) and (29), we shall omit (2), (4), (8), and 
(11), as ineffective in the final equations for d L. (1), (5), (7), and (10) give 

(P) = 1 X 0.8 2 + 1 X l-T + 1 X 1.0 2 + 1 X 0.7 2 = 5.0000, 
and (2), (4), (8), and (11), 

(p) = 1 x 0.03 2 + 1 X 0.01 2 + 1 X 0.08 2 + 1 X 0.01 2 = 0.0075. 

(£) = J_ ' 

(P) 667 ' 

According to (28), we must have, in this case, 

^ 100 

in order to preserve the equivalency of (II.) and (III.), if the original equations (2), 
(4), (8), and (11) are omitted in forming the final equation for d L. We may, then, 
adopt the following system of factors. 

dL. 



No. of 
'iginal Equation. 


Factors in the Method 
of Least Squares. 


Equivalent 
Factors. 


Equivalent 
Factors. 




I. 


II. 


III. 


1 


+ 0.79363 


+ 0.8 


+ 0.8 


2 


— 0.02658 


— 0.03 


0.0 


3 


+ 0.58880 


+ 0.6 


+ 0.6 


4 


+ 0.01318 


+ 0.01 


0.0 


5 


+ 1.73436 


+ 1.7 


+ («- 


6 


— 0.12606 


— 0.1 


— 0.1 


7 


-j- 0.99584 


+ 1.0 


+ 1.0 


8 


— 0.08089 


— 0.08 


0.0 


9 


+ 0.65311 


+ 0.7 


+ 0.7 


10 


+ 0.69957 


+ 0.7 


+ 0.7 


11 


— 0.01315 


— 0.01 


0.0 
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The mean value of H= — for the equations (1), (5), (7), and (10) is 

H = — or < — , 
200 ^ 100 

as required by (28). 

For the second indeterminate dy, the favorable equations are (5), (7), (9), (10). 
The unfavorable ones are (2), (4), (8), and (11). 

(Q) = 10600000, (q) = 9400, 

(?) 1 17 n. ^ * 

— — = — - , ti must be <T — . 

(Q) 1130' ^37 



rf- 



No. of 
Original Equation. 


Factors in the Method 
of Least Squares. 


Equivalent 
Factors. 


Equivalent 
Factors. 




I. 


11. 


III. 


1 


+ 143.66 


+ 0.1 


+ 0.1 


2 


-f- 46.71 


+ 0.05 


0.0 


3 


+ ,358.12 


+ 0.4 


+ 0.4 


4 


-f- 28.39 


+ 0.03 


0.0 


5 


+ 1846.17 


+ 2.0 


— 0.2 


6 


— 227.42 


— 0.2 


7 


+ 1579.03 


+ 2.0 


+ 1.6 


8 


— 67.22 


— 0.07 


0.0 


9 


+ 1329.09 


+ 1.0 


+ (.- 

+ 1.7 


10 


+ 1719.32 


+ 2.0 


11 


— 43.84 


— 0.04 


0.0 







The mean value of H= IB. for the equations (5), (7), (9), (10), is 

Jf=_Lor<I, 

125 ^ 37 

the limit given by (29). 

For d-ir the favorable equations are (1), (3), (5), (9) ; the unfavorable, (2), (4), (11). 

(r) i l 

7-57 = — ? H must be <" — . 

(it) 675 ^92 
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d 7T. 

Equivalent 
Factors. 

III. 
+ 4.0 

0.0 
+ 2.6 

0.0 
— 5.5 
+ 1.0 
+ 0.6 
+ 1.0 
+ 4.0 
+ 1.0 

0.0 
1^ 

92* 

For dcp the favorable equations are (1), (5), 7), (10); the unfavorable, (4), (8), 
(9), (11). 

(£) _J_ 
(S) 465 * 



No. of 
Original Equation. 


Factors in the Method 

of Least Squares. 

I. 


Equivalent 

Facters. 

II. 


1 


+ 0.39493 


+ 4.0 


2 


+ 0.02658 


+ 0.3 


3 


+ 0.26208 


+ 3.0 


4 


— 0.01318 


— 0.1 


5 


— 0.54603 


— 5.0 


6 


+ 0.12606 


+ 1.0 


7 


+ 0.06456 


+ 0.6 


8 


+ 0.08089 


+ 0.8 


9 


+ 0.38994 


+ 4.0 


10 


+ 0.12913 


+ 1.0 


11 


+ 0.01315 


+ 0.1 


The mean value of H 


is — , which exceeds the limit H 

en ' 



The value of H is imaginary, or < — 



i 

J 631 465 ' 

The rejection of the unfavorable equations cannot in this instance be compensated 
for by increasing the accuracy of treatment of the favorable ones. 



d 



<p. 



No. of 


Factors in the Method Equivalent 


Original Equation. 


of Least Squares. Factors. 




i. n. 


1 


+ 0.95920 + 1.0 


2 


— 0.20858 — 0.2 


3 


— 0.85234 — 0.9 


4 


— 0.07861 —0.08 


5 


— 2.05662 — 2.0 


6 


— 0.38939 — 0.4 


7 


+ 1.99545 + 2.0 


8 


— 0.09970 — 0.1 


9 


— 0.08439 — 0.08 


10 


— 1.38787 — 1.0 


11 


+ 0.02929 + 0.03 


The mean value of H 


is — ; it should be as above < 0. 
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27 



Equivalent 
Factors. 

m. 
1 \ ioo7 

— 0.2 

0.0 

— 2.0 

— 0.4 
+ 2.0 

0.0 
0.0 



0.0 
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For dfl the favorable equations are (4) and (11); the unfavorable ones, (7) and (9). 

if). JL 

(T) 344* 

H is again imaginary. In rejecting (7) and (9), we have passed the prescribed limits. 

dn. 



No. of 
Original Equation. 


Factors in the Method 

of Least Squares. 

I. 


Equivalent 

Factors. 

II. 


Equivalent 
Factors. 

m. 


1 


— 0.18856 


— 0.2 


— 0.2 


2 


-f- 0.15946 


4-0.2 


+1 


3 


+ 0.14912 


4-0.1 


+1 


4 


+ 0.91704 


4-0.9 


4- (1—0.08) 


5 


— 0.18833 


— 0.2 


— 0.2 


6 


4-0.17176 


4-0.2 


+ 1 


7 


— 0.06040 


— 0.1 


0.0 


8 


— 0.46359 


— 0.5 


1 

2 


9 


— 0.04305 


— 0.04 


0.0 


10 


+ 0.17130 


+ 0.2 


+1 


11 


+ 1.02138 


4-1.0 


+ (1+0.02) 


ii the favorable equations are (2), (6), (8), and the unfavorable (10). 




(U) 769' 


^ 58 






d, 


i. 




No. of 
Original Equation. 


Factors in the Method 

of Least Squares. 

I. 


Equivalent 

Factors. 

II. 


Equivalent 

Factors. 

III. 


1 


+ 0.17387 


4-0.2 


+ 1 


2 


-f 1.25782 


4-1.0 


+0+0 


3 


4- 0.17775 


+ 0.2 


+5 


4 


4- 0.54365 


4-0.5 


+i 


5 


— 0.17445 


— 0.2 


I 

6 


6 


— 1.35441 


— 1.0 


-o+o 


7 


— 0.33750 


— 0.3 


1 

3 


8 


4- 1.22803 


+ 1.0 


+0+0 


9 


4- 0.34268 


4-0.3 


+ 5 


10 


— 0.08360 


— 0.08 





11 


— 0.27187 


— 0.3 


1 

4 
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The mean value of H is — , the computed limit being H < — . 
Applying these factors to the original equations, we obtain, — 



By the Method of Least Squares. I. 

Coef. of Coef. of Coef. of Coef. of Coef. of Coef. of 
dL. dy. dn. dif. da. di. 

= — 371.09 -}- 5.91569 + 7203.91 — 0.09344 — 2.28516 — 0.34664 — 0.18194 
= — 580104.00 -f 7203.91 + 10834225.00 — 49.06 — 3229.77 — 198.64 — 143.05 
= — 113.45 — 0.09344 — 49.06 -f- 0.71917 + 1.13382 + 0.06400 -f- 0.26341 

= + 268.53— 2.28516— 3229.77+ 1.13382 -f- 12.00340— 0.37137— 0.11762 
= + 94.26 — 0.34664 — 198.64 + 0.06400 — 0.37137 + 2.28215 — 0.36136 

= — 31.81— 0.18194— 143.05+ 0.26341— 0.11762— 0.36136+ 5.62456 

By the Equivalent Factors. II. 

Coef. of Coef. of Coef. of Coef. of Coef. of Coef. of 
dL. dy. dn. dif. da. di. 

= — 372.23 + 5.89945 + 7208.64 — 0.04729 — 2.22332 — 0.33807 — 0.19932 
= — 662.04 + 7.85800 + 11830.88 — 0.20089 — 3.15658 — 0.16440 — 0.48494 
= — 1051.81 — 0.02274 + 126.84 + 7.02511 + 10.51379 + 0.45680 + 3.10941 
= + 137.76 — 1.90276 — 2453.82 + 1.15752 + 11.43969 — 0.33172 — 0.14026 
• = + 84.26 — 0.42389 — 250.17 + 0.05158 — 0.45069 + 2.27706 — 0.39471 
= — 42.98 — 0.20005— 247.16 + 0.29878— 0.04721 — 0.30249 + 4.50959 

By the Equivalent Factors. III. 

Coef. of Coef. of Coef. of Coef. of Coef. of Coef. of 
dL. dy. dn. dip. da. di. 

= — 37L03 + 5.97864 + 7296.25 — 0.06706 — 2.33932 — 0.37876 — 0.08024 
= — 560.48 + 7.21446+10853.11—0.05221— 3.26575 — 0.20615 — 0.18982 
= — 1061.36 — 1.13102 — 959.72 + 7.19887 + 11.91488 + 0.34031 + 3.07535 
= + 251.17 — 2.11274 — 2917.00 + 1.13986 + 11.78434 — 0.37820 + 0.08665 
= + 97.46 — 0.29118— 75.36 + 0.07546— 0.22584 + 2.29772 — 0.39567 
= — 54.37 — 0.12928 — 4.62 + 0.26533 — 0.24548 — 0.37541 + 5.56738 

In order to test the numerical values in these final equations by (33), the decimal 
points for L, M, N, &a, and for the coefficients of d 7, should be changed three 
places to the left. This being done, we have for II. 

2 (P., Qi, Li, M u ) = + 57.87362, 

2 (S 2) — + 57.87354. 

And for III. 

2 (Pi, Qi, £1, M u ) = + 57.56917, 

2 (S 2) = + 57.56916. 
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In both cases, the agreement is as near as could be desired. 

We have, then, the following equations by successive eliminations : — 







Coef. of 


Coef. of Coef. of 




Coef. of 


Coef. of Coef. of 






dL. 


dy. dn. 




d ip. 


do.. di. 


= — 


371.09 + 


5.91569 + 


7203.91 — 0.09344 


— 


2.28516 — 


0.34664— 0.18194 


= —138534.00 


+ 


2458225.00 + 62.13 


— 


510.58 + 


213.84 + 73.45 


= — 


115.81 




4- 0.71612 


+ 


1.11063 — 


0.06392+ 0.25868 


o = + 


25.66 






+ 


9.29213 — 


0.36175— 0.57384 


o = + 


75.23 








+ 


2.22346— 0.37766 


o = + 


17.11 




II. 






+ 5.42383 






Coef. of 


Coef. of Coef. of 




Coef. of 


Coef. of Coef. of 






dL. 


dy. dn. 




dif. 


da. di. 


= — 


372.23 + 


5.89945 + 


7208.64— 0.04729 


— 


2.22332 — 


0.33807— 0.19932 


= — 


166.23 


+ 


2229.06— 0.13790 


— 


0.19504 + 


0.28590— 0.21945 


= — 


1041.71 




+ 7.03449 


+ 


10.51875 4- 


0.43566+ 3.12386 


o = + 


176.07 






+ 


9.01520 — 


0.49449— 0.72094 


o = + 


100.43 








+ 


2.17690— 0.46611 


o- + 


16.00 




III. 






+ 4.24758 






Coef. of 


Coef. of Coef. of 




Coef. of 


Coef. of Coef. of 






dL. 


dy. dn. 




dtp. 


da. di. 


= — 


371.03 + 


5.97864 + 


7296.25— 0.06706 


— 


2.33932 — 


0.37876— 0.08024 


= — 


112.76 


+ 


2048.68 4- 0.02871 


— 


0.44288 4- 


0.25090— 0.09299 


= — 


1108.40 




4- 7.18029 


+ 


11.56325 + 


0.21715+ 3.07926 


= + 


274.45 






+ 


9.07933 — 


0.50448— 0.43778 


o = + 


117.10 








+ 


2.22130— 0.43492 


= + 


30.76 










+ 5.33957 



From which are finally deduced the values of the six unknown quantities : — 



II. 



III. 



di = — 3.15 


— 3.77 


— 5.76 


d a = — 34.37 


— 46.94 


— 53.85 


d q> = — 4.29 


— 22.41 


— 33.50 


d w = + 166.44 


+ 186.17 


+ 212.41 


d y = + 0.054335 


+ 0.08963 


+ 0.05115 


dL = — 3.06 


— 56.32 


— 14.58 



The following are the outstanding errors of the original equations : 
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I. 



II. 



III. 



1 — 111.00 


— 155.51 


— 127.23 


2 — 8.31 


— 3.70 


— 7.29 


3 + 59.18 


-f 59.08 


4- 84.73 


4 — 36.67 


— 47.02 


— 54.55 


5 + 19.92 


+ 21.77 


+ 33.12 


6 + 0.07 


-f 6.92 


4- 19.63 


7 -f 85.77 


-f 54.52 


4- 16.00 


8 + 25.01 


+ 35.37 


4- 38.66 


9 -f 135.88 


-f- 157.61 


4- 144.46 


10 — 204.64 


— 155.64 


— 174.86 


11 + 83.44 


4- 69.63 


4- 64.28 


us of the squares of these 


errors are 




I. 


II. 


III. 


Q = 92919 


88556 


85205 



It is obvious that the solution I., as given by Gauss in the memoir above quoted, 
is incorrect, since the sum of the squares of the errors should be less by the method 
of least squares than by any other mode of combination.* A revised solution gives 
the following equations for I. Those for II. and III. are repeated for the sake of 
comparison. 

I. (Revised solution.) 







Coef. of 


Coef. of Coef. of 


Coef. of 


Coef. of 


Coef. of 






dL. 


dy. dit. 


dif. 


da. 


di. 


= — 


371.09 4- 


5.91569 + 


7203.91 — 0.09344 — 


2.28516 — 


0.34664 — 


0.18194 


o = — : 


128204.00 


+ : 


2061567.00 4-64.72780 — 446.986 4-223.4856 4~ 


78.50985 


= — 


115.28 




4- 0.71566 + 


1.11176 4- 


0.05150 4- 


0.25807 


o = + 


276.46 




+ 


9.29666 — 


0.53681 — 


0.57179 


o = + 


110.68 






+ 


2.20290 — 


0.43212 


= 4- 


41.94 




II. 




+ 


5.40298 






Coef. of 


Coef. of Coef. of 


Coef. of 


Coef. of 


Coef. of 






dL. 


dy. dit. 


dtp. 


da. 


di. 


= — 


372.23 4- 


5.89945 + 


7208.64— 0.04729 — 


2.22332 — 


0.33807 — 


0.19932 


= — 


166.23 


+ 


2229.06— 0.13790 — 


0.19504 + 


0.28590 — 


0.21945 


= — 


1041.71 




4- 7.03449 + 


10.51875 4- 


0.43566 + 


3.12386 


= 4- 


176.07 




+ 


9.01520 — 


0.49449 — 


0.72094 


o = + 


100.43 






+ 


2.17690 — 


0.46611 


= + 


16.00 








+ 


4.24758 



* The errors of equations (10) and (11) given in the Disq. de Elementis Palladis are 
4- 83".01 ; the joint effect of these would increase still further the discrepancy in O. 



-216".54 and 
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III. 



= — 


Coef. of 
dL. 

371.03 + 5.97864 + 


Coef. of 
dy. 

7296.25 — 


Coef. of 
dit. 

0.06706 — 


Coef. of 
d(f. 

2.33932 — 


Coef. of 
dn. 

0.37876 — 


Coef. of 
di. 

0.08024 


= — 


112.76 + 


2048.68 + 


0.02871 — 


0.44288 + 


0.25090 — 


0.09299 


= — 


1108.40 


+ 


7.18029 + 


11.56325 + 


0.21715 + 


3.07926 


= + 


274.45 




+ 


9.07933 — 


0.50448 — 


0.43778 


= + 


117.10 






+ 


2.22130 — 


0.43492 


= + 


30.76 








+ 


5.33957 



Giving, for the unknown quantities, 



(46.) 



I. (Revised 


solution.) 


II. 


III. 


di =— 7.76 




// 

— 3.77 


// 
— 5.76 


d Q = — 51.76 




— 46.94 


— 53.85 


d <p = — 33.20 




— 22.41 


— 33.50 


dn- =+219.19 




+ 186.17 


+ 212.41 


d y = + 0.05401 




+ 0.08963 


— 0.05115 


dL = — 15.68 




— 56.32 


— 14.58 



These results will form the basis for a comparison of the relative probability of the 
three systems. 

In view of the apparently gross discrepancies of the solutions II. and III. from I., 
and with the admitted fact before us that the latter reduces perfectly to a minimum the 
sum of the squares of the residual errors, we should at first sight scarcely be disposed 
to admit that there is in reality no sensible difference in the accuracy of the three 
systems. This unfavorable impression must have its origin in the erroneous infer- 
ence, that, inasmuch as I. is theoretically the best of all solutions, the errors of the 
results of other systems must be proportional to their deviation from I. 

For instance, in I. we have dL = — 16", but in II. dL = — 56"; from which 
is inferred that in II. the error of dL is proportional, or nearly so, to — 16" -f- 56" 
= -j- 40" ; and in a similar way, d 12, dq>, &c. will be condemned or approved by 
the standard of their disagreement or accordance with the corresponding values 
in I. 

Such a mistake may, we repeat, be easily committed, in a cursory glance at the num- 
bers in question. That it is a violation of the fundamental principles of the method 
of least squares will appear on a very little consideration. 

To apply the criterion proposed in that method for determining the relative proba- 
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bility of the three systems, we obtain first the residual errors of the original equa- 
tions as follows : — 

Residual Errors of Original Equations. 

I. (Revised solution.) II. III. 

1 — 125.45 — 155.51 — 127.23 

2 — 9.14 — 3.70 — 7.29 

3 -f- 86.62 -f 59.08 -f- 84.73 

4 — 53.78 — 47.02 — 54.55 

5 + 32.12 + 21.77 + 33.12 

6 + 22.92 + 6.92 -f 19.63 

7 -f 21.00 -f 54.52 -f 16.00 

8 + 35.51 -f 35.37 -f 38.66 

9 + 149.37 -f- 157.61 -f- 144.46 

10 — 169.77 — 155.64 — 174.86 

11 + 66.95 + 69.63 -f 64.28 

The sums of the squares of these errors are 

I. II. III. 

Q = 85091 88556 85205 

The probable value of a residual error for one of the original equations — since there 
are eleven equations* and six indeterminates — will be obtained from the expression 

H= 0.67459 J — — , 
Nfll — 6 

giving the values 

I. II. III. 

(47.) /* = ± 88".00 ± 89".77 ± 88".06. 

The agreement of /t in I. and III. is a conclusive proof of the equivalency of the two 
solutions, notwithstanding the freedom which has been exercised in the choice of fac- 
tors for the latter. In the case of II., the discrepancy amounts to - of /* ; but it will 
be noticed that the factors actually employed were based upon a system of representa- 
tion considerably less exact than the series (2). It might easily be shown that this 
difference would have been reduced to less than one fifth of its present amount if the 
numbers for the factors had been selected from (2), agreeably to the conditions by 
which II. has been denned. 

* Equation 9 having been excluded from each of the solutions. 
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If we compute the limits of these values of n, we find that it is only an even 
chance that /*, for the best solution, is comprised within the limits 

p= ± 69".3, and /» = ± 106".7. 

Such is the extreme uncertainty of the only element by which the question of pref- 
erence between I., II., and III. can be decided.. Compared with it, the inconsiderable 
differences which we find between the values of //. in I., II., and III. will admit of 
but the single inference, that either of the systems of values of di, dfl, dq>, &c, 
presented in (46), notwithstanding their great disagreements with each other, actu- 
ally fulfils the criterion of accuracy proposed in the method of least squares so 
nearly, that it is impossible to give a decisive reason for adopting one rather than 
another as the most probable solution. It is worthy of notice, moreover, that the 
arithmetical mean of the above residual errors, irrespective of their signs, is less 
in II. and III. than in I. Thus in this instance the latter would rank lowest of 
the three, if we were to compute the relative probabilities according to a process 
recommended by the highest authorities* as the most suitable for ordinary use, in 
which the probable error is directly proportional to the arithmetical mean of the errors 
irrespective of their signs. 

* Laplace, Theorie Analytique des Probdbilitis ; Gauss, Zeitschrift fiir Astronomie, B. I. ; Peters, Astr. 
Nach., No. 1034. 



